Multi-State VALBOND for Atomistic Simulations of Hypervalent Molecules, Metal Complexes, and Reactions.

The implementation, validation, and application of the multi-state VALBOND method for transition-metal-containing and hypervalent molecules are presented. This approach is particularly suited for molecules with unusual shapes and systems that need to be described by a superposition of resonance structures, each of which satisfies the octet rule. The implementation is based on the original VALBOND force field and allows us to smoothly switch between resonance structures, each of which can be characterized by its own force field, including varying charge distributions and coupling terms between the states. The implementation conserves total energy for simulations in the gas phase and in solution and is applied to a number of topical systems. For the small hypervalent molecule ClF3, the barrier for pseudorotation is found to be 4.3 kcal/mol, which compares favorably with the experimentally measured value of 4.8 kcal/mol. A transition-metal-containing complex, cisplatin, is characterized by six resonance states, for which the vibrational spectrum is found to be in good agreement with experiment. Finally, umbrella sampling simulations of the SN2 reaction BrMe + Cl- → Br- + MeCl in solution yield a barrier height of 24.6 kcal/mol, in good agreement with experiment (24.7 kcal/mol).

A dual anchoring strategy for the localization and activation of artificial metalloenzymes based on the biotin-streptavidin technology.

Artificial metalloenzymes result from anchoring an active catalyst within a protein environment. Toward this goal, various localization strategies have been pursued: covalent, supramolecular, or dative anchoring. Herein we show that introduction of a suitably positioned histidine residue contributes to firmly anchor, via a dative bond, a biotinylated rhodium piano stool complex within streptavidin. The in silico design of the artificial metalloenzyme was confirmed by X-ray crystallography. The resulting artificial metalloenzyme displays significantly improved catalytic performance, both in terms of activity and selectivity in the transfer hydrogenation of imines. Depending on the position of the histidine residue, both enantiomers of the salsolidine product can be obtained.

Toward a Broadly Applicable Force Field for d(6)-Piano Stool Complexes.

Three-legged piano stool complexes are prototypical organometallic complexes relevant to a wide range of chemically relevant questions. Force field parametrization of transition-metal complexes is difficult and underdeveloped, and metal-specific force fields and software are required. Here we report our efforts to derive parameters for the conventional CHARMM and the Valbond-CHARMM force fields for d(6)-piano stool complexes. In Valbond-CHARMM, the usual angular term is replaced with hybrid orbital strength functions. These functions describe the energy not only of distorted bond angles around the minimum but also at very large distortions. Structure optimizations led to a good agreement between the calculated force field and the X-ray structures. They were comparable to RMSDs obtained between X-ray and DFT structures. In addition, and contrary to treating the systems with DFT, molecular dynamics simulations on the multiple nanosecond time scale are possible and allow to compute meaningful structural and energetic observables. Explicit solvent simulations of the complexes in methanol and water allow to determine the solvent distribution around the complexes. The parametrization presented here will be a useful starting point for dynamics investigations of catalysts in structurally more demanding environments.

Identification of two-histidines one-carboxylate binding motifs in proteins amenable to facial coordination to metals.

Among natural metalloenzymes, the facial two-histidines one-carboxylate binding motif (FTM) is a widely represented first coordination sphere motif present in the active site of a variety of metalloenzymes. A PDB search revealed a total of 1685 structures bearing such FTMs bound to a metal. Sixty statistically representative FTMs were selected and used as template for the identification of structurally characterized proteins bearing these three amino acids in a propitious environment for binding to a transition metal. This geometrical superposition search, carried out using the STAMPS software, returned 2320 hits. While most consisted of either apo-FTMs or bore strong sequence homology to known FTMs, seven such structures lying within a cavity were identified as novel and viable scaffolds for the creation of artificial metalloenzymes bearing an FTM.

Molecular Mechanics Force Field for Octahedral Organometallic Compounds with Inclusion of the Trans Influence.

Efficient calculation of the properties of metal-containing complexes relevant to catalysis is of major interest for better characterizing and optimizing the catalysts. For this, a new force field, called VALBOND-TRANS here, is proposed. It is based on the existing VALBOND force field of Landis and co-workers, extended by adding terms that account for electronic effects such as the trans influence of ligands on bond lengths and relative energies. Parameters and results for model octahedral complexes of Ru, Os, Rh, and Ir are determined and discussed. The model is then applied to the study of reactive intermediates involved in asymmetric hydrogenation catalyzed by iridium complexes with chiral phosphinooxazolines (PHOX) ligands. The new force field explores and capitalizes on the separation of electronic and steric effects on the stability of different diastereomers and reproduces DFT results which are consistent with experimental observations.